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Abstract 

A recent work on the resummation of fermionic in-medium ladder diagrams to all or- 
■ ders is extended by considering the effective range correction in the s-wave interaction 

f-*) . and a (spin-independent) p-wave contact-interaction. A two-component recursion gen- 

CN \ erates the in-medium T-matrix at any order when off-shell terms spoil the factorization 

of multi-loop diagrams. The resummation to all orders is achieved in the form of a ge- 
ometrical series for the particle-particle ladders, and through an arctangent-function for 
the combined particle-particle and hole-hole ladders. One finds that the effective range 
correction changes the results in the limit of large scattering length considerably, with the 
effect that the Bertsch parameter ^„ nearly doubles. Applications to the equation of state 
of neutron matter at low density are also discussed. For the p-wave contact-interaction 
the resummation to all orders is facilitated by decomposing tensorial loop-integrals with a 
^ ■ transversal and a longitudinal projector. The enhanced attraction provided by the p-wave 

. ladder series has its origin mainly in the coherent sum of Hartree and Fock contributions. 

T— I ■ 

>■ PACS: 05.30.Fk, 12.20.Ds, 21.65+f, 25.10. Cn 
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1 Introduction and summary 

O \ Dilute degenerate many-fermion systems with large scattering lengths are of interest, e.g. for 
^ \ modeling the low-density behavior of nuclear and neutron star matter. Because of the possibility 
to tune magnetically atomic interactions, ultracold fermionic gases provide an exceptionally 
valuable tool to explore the non-perturbative many-body dynamics involved in the crossover 
^ . from the superconducting to the Bose-Einstein condensed state (for a recent comprehensive 
^! review of this fascinating field, see ref. [1]). Of particular interest in this context is the so- 
called unitary limit in which the two-body interaction has the critical strength to support a 
bound-state at zero energy. As a consequence of the diverging scattering length, a — )■ oo, 
the strongly interacting many-fermion system becomes scale- invariant. Its ground-state energy 
is then determined by a single universal number, the so-called Bertsch parameter ^, which 
measures the ratio of the energy per particle E(kf)^°°^ to that of a free Fermi gas E(kf)^^'> = 
Skj/IOM. Here, kf denotes the Fermi momentum and M the fermion mass. 

The calculation of ^ is an intrinsically non-perturbative problem which has been approached 
by numerical quantum Monte Carlo simulations. As the state of the art, a value of ^ ~ 0.38 
emerges at present [2, 3]. It is consistent with a recent precision measurement of the equation 
of state by the MIT group, which gives ^ = 0.37 ±0.01 [4] (extrapolating to zero temperature). 
Among the analytical approaches let us mention the e-expansion by Nishida and Son [5], where 
e = 4 — d with d the number of space dimensions. A naive extrapolation to d = 3 gives 
^ = 0.475, but with inclusion of higher order corrections and Fade interpolations between the 
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(otherwise diverging) expansions around d = 2 and d = 4, a. good value of ^ = 0.38 ±0.01 could 
be obtained (see chapter 7 in ref.[l]). Because of the very large neutron- neutron scattering 
length ann — 19 fm, neutron matter at low densities is supposed to be a Fermi gas close to 
the unitary limit. Quantum Monte Carlo simulation [6, 7] and other sophisticated many-body 
calculations of neutron matter give indications for a value of — 0.5. 

In a recent work [8] the complete resummation of the (combined) particle-particle and hole- 
hole ladder diagrams generated by a contact-interaction proportional to the scattering length 
a has been achieved. A key to the solution of this (restricted) problem has been a different 
organization of the many-body calculation from the start. Instead of treating (propagating) 
particles and holes separately, these are kept together and the difference to the propagation in 
vacuum is measured by a "medium- insertion" . The following identical rewriting of the (non- 
relativistic) particle-hole propagator: 



explains the basic principle of the approach. In that organizational scheme the pertinent in- 
medium loop is complex- valued, and therefore the contribution to the energy per particle E{kf ) 
at order a" is not directly obtained from the (n — l)-th power of the in-medium loop (see Fig. 1). 
However, after reinstalling the symmetry factors l/(j + 1) which belong to diagrams with j + 1 
double medium-insertions, a real- valued expression is obtained for all orders a". Known results 
about the low-density expansion [9, 10] up to and including order could be reproduced with 
improved numerical accuracy. The emerging series in akf could even be summed to all orders 
in the form of a double-integral over an arctangent-function. In that explicit representation the 
unitary limit, a oo, could be taken straightforwardly and the value ^ = 0.507 has been found 
for the Bertsch parameter [8]. This number is to be compared with the value ^^^^^ ~ 0.237 
resulting from a resummation of particle-particle ladder diagrams only [1 1] . Clearly, since only 
the (infinite) subset of particle-hole ladder diagrams [8] is computed in the normal phase one 
cannot expect accurate results for the unitary Fermi gas, whose true ground state is a superfluid. 
But interestingly, extrapolations of the equation of state of the unitary Fermi gas from finite 
to zero temperature, which smoothly pass over the pairing transition at give indications for 
a Bertsch parameter — 0.5 in the normal phase [1, 4, 12]. In comparison to this synthetic 
result, the value = 0.507 of ref.[8] appears to be a rather good finding. 

The purpose of the present paper is to extend the resummation method of ref.[8] by in- 
cluding the effective range correction in the s-wave interaction and by considering also a 
spin-independent p-wave contact interaction to all orders. The effective range correction to 
a contact-interaction (proportional to the scattering length a) has been studied first by Schafer 
et al. in ref.[ll] for the resummed particle-particle ladder series. The imphcations of this 
correction on the equation of state have, however, not been explored in much detail and only 
rough perturbative estimates have been given. In a related work by Schwenk and Pethick [13] 
it has been argued that the inclusion of the effective range = 2.75 fm leads to an improved 
description of neutron matter at low densities. Their implementation of the effective range rg 
into the particle-particle ladder series assumes a form similar to the effective range expansion 
for neutron- neutron scattering in free space. However, a proper diagrammatic calculation [11] 
shows that in the medium the combination of scattering length a and effective range pro- 
duces additional (density-dependent) off-shell terms which spoil the factorization of multi-loop 
diagrams. A resummation of the particle-particle ladder diagrams to all orders in the form of a 
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geometrical series is still possible, but it differs from the expression suggested by the effective 
range expansion in vacuum. 

The present paper is organized as follows: In section 2 we perform the resummation of 
fermionic in-medium ladder diagrams for a contact-interaction that includes the effective range 
Tg. It is demonstrated that a two-component recursion generates the in-medium T-matrix 
at any order when off-shell terms spoils the direct factorization of multi-loop diagrams. The 
resummation to all orders is still achieved in the form of a geometrical series for the particle- 
particle ladders, and through an arctangent-function for the combined particle-particle and 
hole-hole ladders. The second more comphcated case involves yet a (combinatorial) conjecture 
because the arctangent-series could be verified so far with diagrammatic methods only to some 
finite order. One finds that the effective range correction changes the results in the limit 
of large scattering length a ^ oo considerably, with the effect that the Bertsch parameter 

(for the normal phase) increases. For the combined particle-hole ladder series one gets 
ini^s 7^ 0) = 0.876, and for the resummed particle-particle ladder series alone the value is 
^n^K^'s 7^ 0) ~ 0.43. Applications of these resummation results to the equation of state of 
neutron matter at low density are presented and discussed in comparison. Section 3 deals with 
the analogous resummation to all orders for a (spin-independent) p-wave contact-interaction. 
In this case the factorization of multi-loop diagrams is achieved more directly by decomposing 
tensorial loop- integrals into a transversal and a longitudinal part. These two parts decouple and 
can therefore be resummed separately through an arctangent-function or a geometrical series. 
The enhanced attraction provided by the p-wave ladder series (in the limit of large scattering 
volume af oo) has its origin mainly in the coherent sum of Hartree and Fock contributions. 
In the appendix analytical results for the resummation of particle-hole diagrams to all orders 
are presented for a two-component fermionic system (such as nuclear matter) with two different 
scattering lengths, and at, in the presence of an isospin-asymmetry, kn > kp. A special case 
of such an asymmetric system is the (unitary) Fermi gas in the spin-unbalanced configuration. 

Presumably the results of the present paper arc contained (to a certain extent) as special 
cases in the work of Lacour et al. [14]. In that paper non-perturbative methods for chiral 
effective field theory of finite density nuclear systems have been developed. The approach is 
based on a partial wave representation of leading order nucleon-nucleon interactions which are 
summed to all orders in the medium. A detailed comparison of the present results (obtained 
from simple contact-interactions) with those of ref.[14] is certainly interesting, but this requires 
an extensive communication between both groups in order to translate the different formalisms 
correctly into each other and to establish the common range of applicability. 

The study of strongly-coupled fermionic many-body systems (such as neutron matter) at low 
density with quantum Monte Carlo simulations [15, 16, 17] is presently an active research area. 
The analytical results of the present work, which of course cover only a restricted part of the 
many-body dynamics, could be potentially useful for interpreting these numerical simulations. 

2 Resummation with inclusion of the effective range 

In this section we develop methods which allow for a complete resummation of fermionic in- 
medium ladder diagrams for an s-wave contact-interaction that includes the effective range 
correction. We recall from rcf.[8] that after opening a closed multi-loop diagram (represent- 
ing the energy density) at adjacent medium- insertions a planar ladder diagrams results (see 
Fig. 1). If the contact-interaction is momentum- independent all these (open) multi-loop dia- 
grams factorize into powers of the in- medium loop. In the ordering scheme implied by eq.(l) 



3 



<^AAAAA• 
<^AAAAA• 
<^AAAAA<. 

"fWW\A» 

'\AAAAAii 
'^AAAA/^l 

Figure 1: Left: Closed multi-loop diagram representing a contribution to the energy density. 
Right: Planar ladder diagram obtained by opening at a minimal pair of adjacent medium- 
insertions. Wiggly lines symbolize contact-interactions. 
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Figure 2: In-medium loop organized in the number of medium-insertions (=). The central 
diagram has a reflected partner. The external momenta |pi^2| < ^/ are inside the Fermi sphere. 
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the in-medium loop is composed of contributions with zero, one, and two medium-insertions 
(see Fig. 2). The first part represents the rescattering in the vacuum and the other two parts 
incorporate systematically the (Pauh-blo eking) effects due to the fermionic medium. The re- 
sulting in-medium loop is complex- valued, but after a careful combinatoric analysis of the 
closed multi-loop diagrams a real- valued expression for the energy density is obtained at each 
order. In the next step the complete rcsummation to all orders can be performed through an 
arctangent-function that depends on both the real and imaginary part of the in-medium loop. 
The aim is to generalize this novel approach to an s-wave contact- vertex: 

Co + C^iql + Cd. (2) 

which includes the effective range correction. Here, gin and g^ut denote half of the momentum 
difference in the initial and the final state, respectively. Such a form of the contact-vertex is 
actually dictated by Gahlei invariance. The relations of the couphng constants Cq and C2 to 
the scattering length a and the effective range parameter are: 

Note that wc arc choosing the sign-convention such that a positive scattering length a > 
corresponds to attraction. The repeated rescatterings via the contact- vertex cq.(2) in vacuum 
(i.e. the left diagram in Fig. 2 iterated to all orders) leads to an s-wave phase shift 5o(g) of the 
form: 

2 

taiiSo{q)=aq (4) 

with q the ccnter-of-mass momentum. We are using dimensional regularization, where all (scale- 
free) power divergences are set to zero: dll"' = 0. The real part of the in-medium loop is 
generated by the central diagram in Fig. 2. The pertinent integration formula for the real part 
reads: 

d^l -M 
(2^ 1 2 - g-: 

with the factors from the lower and upper interaction-vertex to be inserted at the place of the 
dots. Here, we have introduced the half sum P = [pi +p2)/'2 and half difference q = (pi — p2)/2 
of the external momenta |pi_2| < ^/ inside the Fermi sphere. All three diagrams in Fig. 2 
contribute to imaginary part of the in-medium loop and the pertinent integration formula 
reads: 

/(0^ ~ ^^^^ " " ^"^'^ ^^^^ - |p + ri) . . . . (6) 

Note that we have left out a term [1 — 9{. . .)][! — 9{. . .)] which vanishes on-shell (/^ = q'^) 
due to energy conservation and Pauli-blocking (sec cq.(8) in ref. [8]). Let us now analyze more 
closely the one-loop diagrams in Fig. 2 with two momentum- dependent Co,2 vertices. For the 

— * — * 

lower vertex the assignment q\a = q, (fcut = I with I the loop-momentum holds, while for the 
upper vertex it is gin = ^, ?out — q- Actually, if a further loop is considered the second g^ut is set 
to the new loop momentum / ' and therefore this second g^ut should not be fixed immediately. 
One recognizes that the product of interaction vertices contains terms proportional to and 
l^. With inclusion of the energy denominator in eq.(5) these can be decomposed into on-shell 
vertices (/^ — )■ g*^) and polynomial off-shell pieces: 

^-^ = J^ + 1, ^^ = ^^ + P + g-. (7) 
l'^ — q^ l'^ — (p P — q^ l"^ — q^ 
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These off-shell terms produce additional real-valued contributions which enter into the next 
loop, where by the same mechanism new off-shell terms are generated, and so on. The fac- 
torization of multi-loop diagrams into powers of the in-medium loop gets thus spoiled by the 
presence of the C2 coupling proportional to the effective range r^. Although it looks compli- 
cated at first sight, one can recursively follow the continuous proliferation of off-shell terms 
by simultaneously keeping track of two quantities: the (on-shell) in-medium T-matrix T„ and 
the coefficient of g^^^. Injecting the expression T„ -|- (3n{(lout ~ order n with g^ut ^ i 

and concatenating it through the in-medium loop with one additional Co,2 vertex produces as 
an output the same expression at order n + 1. As a transparent result one deduces a linear 
two-component recursion of the form: 

(8) 

with the starting values Ti = Cq + 2C2q^ and (3i — C2. The pertinent 2x2 matrix A4 reads: 

. . _ /(Co + 2C2?2)£ + C2X (Co + C2q^)X + C2Y\ 
•^-[ C2L C2X ) ' 

with X, Y the off-shell terms originating from Fermi sphere integrals over the terms 1, P in 
eq.(7) and C the complex-valued in-medium loop (see eqs. (14-17)). Although the polynomial 
expressions for the in-medium T-matrices T„ get increasingly more complicated with increasing 
n, their total sum can be easily calculated with the help of matrix-inversion: 

E Tn = (1, 0) [1 - m] ^^^'^'^ = [Kff^ = E V:, C--' . (10) 

The second equality reveals that the reciprocal of this infinite scries splits off instantaneously 
the complex- valued in-medium loop C and the third equality shows that the result is identical 
to a geometrical series with an effective density and momentum dependent potential: 

Co + 2C2g^ + C|(y - Xg^) 

(1 - C2xy ■ ^^^^ 

The beginning and the end of eq.(lO) represent a remarkable identity. It means that although 
the individual multi-loop diagrams did not factorize, their total sum agrees exactly with a 
geometrical series one would have obtained by using the effective potential V^s and assuming 
factorization. Note that V^s depends in a non-linear way on the off-shell term X and when 
dropping the effective range C2 = one recovers \4ff = Co = Ana/M. A resummation of the 
in-medium T-matrix with inchision of all partial waves has been proposed in section 4.2 of 
ref.[14] (assuming on-shcll factorization). 

The intermediate resummation result in eq.(lO) does not yet provide the correct integrand 
for computing the energy per particle E{kf) because it is complex- valued. Following the in- 
structive diagrammatic analysis in section 4 of ref.[8] we make the conjecture that the proper 
resummation of in-medium ladder diagrams to all orders is achieved by the substitution: 

1 1 Im£ 

T :7 — >■ T — ;7 arctan — — - . 12 

We have verified by a detailed diagrammatic analysis that the arctangent-prescription in eq.(12) 
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Figure 3: In- medium ladder diagrams contributing to the energy density at fourth order. The 
third and fifth diagram have symmetry factors 1/2 and 1/4, respectively. 

leads to the correct result in a perturbative expansion in the coupling constants Co,2 (at least) 
up to sixth order. In order to present an example, the corresponding set of diagrams at fourth 
order is shown in Fig. 3. These decorated diagrams have symmetry factors 1, 1, 1/2, 1 and 1/4, 
in the order shown. Taking these symmetry factors into account the integrand for computing 
the energy density at fourth order has the form: 

Tl+T;Ti{2ilmC) + ^{T;)\2ilmC) + T;{Ti)\2ilmCy + ^{Ti)\2ilmCf , (13) 

where T* denotes the complex-conjugate^ of the n-th order T- matrix T!„ to be calculated 
recursively via eq.(8). This sum of five complex terms is indeed real- valued and it agrees 
exactly with the expression one obtains by expanding the arctangent- function in eq.(12) to 
fourth power in the couplings Co,2- The relevant real- valued integrands at second and third 
order are: T* + |(Ti)2(2iIm£) and T* + T*Ti{2ilmC) + ^{Tif{2ilmCy. Going step by step 
perturbatively to higher n is in principle straightforward, it requires to check in each case the 
equality of two very long polynomial expressions. A complete verification of eq.(12) to all orders 
n requires the consideration of all possible partitions (^i, n2, . . . ,nj) of n = rii + n2 + . . . + fij 
together with their symmetry factors. The sequence (ni, n2, . ■ ■ ,nj) counts here the numbers 
of contact-interactions between consecutive double medium- insertions (see Fig. 3). Since the 
partitions of n are in one-to-one correspondence to the inequivalent irreducible representations 
of the permutation group Sn there might exist some sophisticated combinatorial identities, 
which could help to solve this problem completely. 

Before turning to the final result for the resummed energy per particle we give the explicit 

^We are using here the fact derived in section 3 of rcf.[8] that the complex- valued in- medium loop C goes 
over into its complex-conjugate if the diagram with double medium insertions is taken out. 
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analytical expressions for the off-shell terms: 



^ = -1^. = -1^(5^^ + 3). (14) 



and for the complex-valued in-medium loop: 

. Mk 



f { +i7r7(s,K)}. (15) 



47r2 

The dimensionless variables s = \pi + p2\/2kf and k. = \pi — p2\/2kf fulfil the constraint 
s'^ + k'^ < 1, since |pi,2| < kf. The dimensionless real part R{s,k) and imaginary part I{s,k) 
introduced in eq.(15) have the form [8]: 

H(s. «) ^ 2 + 1|1 - (. + .f] In + - (- 1- 1^ . (16) 

I{s, k) = «; ^(1 - s - k) + ^(1 - - K^) e{s + K-l). (17) 

The same expression for the in- medium loop has been derived in appendix C of ref.[14]. With 
these ingredients the resummed energy per particle is readily calculated. The 2-ring Hartree 
and 1-ring Fock diagram add up with their spin- factors as 4 — 2 = 2. One uses the reduction 
formula cq.(15) in ref.[8] for the integral over two Fermi spheres \pi,2\ < kf to cancel yet the 
factor 1/Im£ in eq.(12), and obtains the following expression for the interaction energy per 
particle: 







with the auxiliary function: 



= (raSfplTP - ) + ^(^^ - - • 

depending on the physical parameters: scattering length a, effective range Vs and Fermi momen- 
tum kf. When dropping the effective range one recovers Q{a, 0) = akf/i: and in the limit of large 
scattering length, lini f2(a, r^) = 3(5/«^—5s^—3)/20 becomes independent of r^A;/. This compar- 
ison demonstrates that the limits a — > 00 and — > do not commute. The arctangent-function 
occurring in eq.(18) refers to the usual branch with odd parity, arctan(— x) = — arctanx, and 
values in the interval [— 7r/2, 7r/2]. We note that the leading order correction involving the 
effective range, agrees with the low-density expansion of ref.[9]. 

Since the unitary limit a — )■ 00 can be performed straightforwardly for eq.(19) a first interest 
aims at the Bertsch parameter. A numerical evaluation of the double-integral in eq.(18) gives 
Cn{i^s 7^ 0) = 0.876, which is considerably larger than the value ^„(rs = 0) = 0.507 for vanishing 
effective range. The reason for this feature are the off-shell terms X,Y in eq.(ll) which change 
the reciprocal effective interaction V^^^ drastically. Note that the ratio C2/C0 = —arg/^ is 
large. The resummation formula given in eqs.(18,19) offers of course the possibility to study 
various other strong coupling limits a — > 00, by imposing a certain relation between a and r^. 
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2.1 Particle-particle ladders only 

In this subsection we adapt the previous methods to the particle-particle ladders series with 
inclusion of the effective range correction [11]. The particle-particle ladder diagrams are re- 
summed via a geometrical series {V^^—C)~p, in which the quantities X, Y and C take on a differ- 
ent form. For particle-particle intermediate states only the integration region in eq. (5) changes 

— * — * 

effectively from the sum to the union of two partly overlapping Fermi spheres |P ± Z | < kf. 
Performing these integrals, one finds for the pp off-shell terms: 

Mk"^ Mk^ 
^pp = -1^(1 + - s) , = -^(1 + ^)'(6 - 3s + s^) , (20) 

in the region < s < 1 of interest, and for the pp in-medium loop: 

^pp ~ ~ ^pp(^-> 1 l^-^) 
with the familiar particle-particle bubble function (s^ -\- <1): 

^ , N ^ l-\-S-\-K 1 n o, , (1 S)^ - 

1 + S — K 2s 1 — — K'' 

Putting these pieces together, one derives from the resummed particle-particle ladder diagrams 
the following expression for the energy per particle: 



Elkf)^-"^ [dss^ I dKK—- i^'"''' , r|>, (23) 

^ ^ ' Ml I Q (n r 1-1 -I- F (^o i^W ' ^ ^ 



M J npp{a,rs)-'^ + Fpp{s,K) 



with the auxiliary function: 

3ak 



^-^"'^^^ = [12. - a-r,k%l + sn2 - s)]- {''"^' " "^^^^'^ 



Wr^fil + s) 



(2 -5)^2 + 1^(35- 6- s^) 



(24) 



depending on the physical parameters a,rs and kf. Again, when dropping the effective range, 
flpp{a, 0) = akf/n, and in the strong coupling limit Qpp{oo.rs) = 3[(1 + s)(3s — 6 — s^) -|- 5(2 — 
s)fi;2]/[5(l + s)'^(2 — s)2] becomes independent of r^/c/. Because of zeros in the denominator, 
the double- integral in eq.(23) is to be understood as a principal- value integral. In the actual 
numerical treatment we work with a cutting- function, Cut(x) = a; for < A^, Cut(a;) = for 

> A^, in order to implement the symmetrical excision around the first-order poles and study 
the convergence behavior for large A^. It is important to note that our result agrees exactly 
with eq.(13) in ref.[ll] (after matching the conventions for the coupling constants 6*0,2)- The 
underlying structure of a geometrical series with an effective potential V^^^ and the reduction 
to an integral over the quarter unit-disc has seemingly not been recognized in that paper. Note 
that ref.[ll] has investigated in addition the summable series of the particle-hole ring diagrams. 
For this subclass of diagrams (as well as for the isolated holc-holc ladder series) the energy per 
particle E{kf) diverges in the strong coupling limit a ^ 00 (see sections 3 and 4 in ref.[ll]). 

In the unitary limit a — > 00, a numerical integration gives for the Bertsch parameter 
^^PP\rs 7^ 0) ~ 0.43 which amounts to almost twice the value $,^\rs = 0) ~ 0.237 at van- 
ishing effective range. The approach of Schwcnk and Pcthick [13] has employed the form 
Opp(a, rs)~^ — TT{l/akf + rskfK^/2) as suggested by the effective range expansion in free space. 
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neutron density [fm ] 

Figure 4: Energy per particle of neutron matter as a function of the neutron density p„ = 
/cy/STT^. The line labelled APR stems from the sophisticated many-body calculation of ref.[18]. 

2.2 Application to low-density neutron matter 

As an application we consider the equation of state of neutron matter. Due to the very large 
neutron-neutron scattering length a„„ = (18.95 ± 0.40) fm [19] neutron matter at low densities 
is supposed to be a Fermi gas close to the unitary limit. Fig. 4 shows the energy per particle 
of neutron matter as a function of the neutron density pn = kj/Sn^. The full and dashed line 
result from the resummation of (combined particle- hole) ladder diagrams eq.(18) with inclusion 
of the empirical effective range parameter = (2.75 d= 0.11) fm and by dropping it, = 0. 
The resummation of the particle-particle ladders only leads to dashed-dotted line in Fig. 4 and 
the result of the approach by Schwenk and Pcthick [13] is additionally included. In each case 
the sum of kinetic and interaction energy per particle 3A;^/10M -|- E{kf) is plotted. The curve 
labeled APR stems from the sophisticated many-body calculation by the Urbana group [18], 
to be considered as a representative of realistic neutron matter calculations. In fact a recent 
calculation of neutron matter in chiral effective field theory with inclusion of subleading three- 
nucleon forces etc. leads to a very similar equation of state [20] . Only the lower two lines in 
Fig. 4 give a reasonable reproduction of the APR curve at low densities Pn < 0.2 fm~^. Since 
the dimensionless variable Unnkf reaches values up to 34.3 in this density region the behavior 
of the total energy per particle Skj/IOM + E{kf) is almost entirely determined by the value of 
the Bertsch parameter Note that for a correct reproduction of the neutron matter equation 
of state at low densities it has to be close to ^„ ~ 0.5. The additional repulsive effects arising 
from hole-hole and mixed particle-hole ladders are recognizable as the difference between the 
full and dashed-dotted line in Fig. 4. Despite the deviations from the APR curve visible in 
Fig. 4, we reproduce the results of Schwenk and Pethick for densities p„ < 0.02 fm~^ actually 
considered in ref.[13]. 



10 



3 Resummation of p-wave contact interaction 



In this section we extend the resummation method to a spin-independent p-wave contact- 
interaction.^ The pertinent contact-vertex obeying Gahlei invariance reads: 



(25) 



with the couphng constant Ci related to the p-wave scattering volume af by Ci = 127raf/M. 
The repeated rescatterings via this contact-vertex in the vacuum (i.e. the left diagram in Fig. 2 
iterated to all orders) leads to a p-wave phase shift of the form: 



tan5i(5) = {aiqf , 



(26) 



with q the center-of-mass momentum and we have used the rule dW^ — of dimensional 
regularization. In the in-medium loop two consecutive Ci vertices introduce the factor (g^n • 

— * — * 

I) {I- Qout)- We decompose the emerging tensorial integral over lilj into the transversal projector 
Sij — PiPj and the longitudinal projector PiPj, where P = [pi + p2)/2 (see Fig. 2). Since these 
are two orthogonal projectors the factorization of multi-loop diagrams holds separately for the 
transversal and longitudinal part of the tensor to be contracted with ql^qi^^. In this decoupling 
procedure one encounters a transversal in-medium loop whose real and imaginary part read: 



- - - + 2«2 + ^ ' + -f- ^ + «2 - Ask - 1 

3 4 As^ 16s3 



, l+S+K 
X In-r^ r + 



[s - nf - If 



\1 — s — k\ 



16s3 



(s^ + ^ Ask - 1) In 



1 + s 



K 



1 — S -\- K 



I±{s, k) ^ k^ 9{1 - s - k) + 



1-s^-k' 
16s3 



12s^k' - (1 - - Ky 9{s + K-l), 



(27) 
(28) 



as well as a longitudinal in-medium loop whose real and imaginary part read: 



R\\ (s, k) 



3 

X In 



+ 2k^ 



[1- K 



2\2 



2s2 



+ 



K^ 



8s3 



S — K\ 



+ 



[l-S^- K 

8s3 



2\3 



+ k' 



111- 



K 



S + K 



(29) 



I\\{s,k) = k^ 9{1 - s - k) + 



2\3 



[S + K-I). 



(30) 



The diagrammatic analysis in section 4 of ref.[8] emphasizing the important role of symmetry 
factors instructs that for the energy density the proper resummation to all orders is achieved 
through an arctangent-function. By using the operator identity for projectors, arctan(nxXx + 
IIiixii) = n^arctanxj, -|- IIn arctanx||, the transversal and longitudinal parts get separated 
completely. A special feature of the (resummed) p-wave contact-interaction is that the 2- 
ring Hartree and 1-ring Fock diagram add up with their spin-factors as 4 + 2 = 6. While 
gin = (Pi ~ P2)/2 = q for both, there is a sign difference in g^ut = ig* for the Hartree and 
Fock contributions. Putting all pieces together one finds from the resummed ladder diagrams 



^Additional spin-orbit and tensor forces are necessary to describe the splitting of the three "^-Po,i,2 channels. 
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with a (spin-independent) p-wave contact-interaction the following expression for the energy 
per particle: 



E{kf) = -y [dss"^ [ (i«;«;|2arctan7 — , v ■f'^^^'^^r. / 7 

-|- arctan — , , „ — — r > . (31) 

(ai%)-3 + 7r-ii?||(s,«;) J ^ ^ 

The factor 2 is reminiscent of two transversal versus one longitudinal degree of freedom. Another 
interesting feature is that the imaginary parts of the transversal and longitudinal loop functions, 
/^(s, k) and /||(s, k) in eqs.(28,30), appear (with a prefactor s^k) as weighting functions in the 
reduction of integrals over two Fermi spheres |pi,2| < 

As a first application let us consider the low-density expansion of E{kf) in eq.(31), which 
rG3jCis * 

F,/, X r 9 , , ,o 1033- 156 ln2, , 1 

The linear term in the scattering volume agrees with ref.[9] (after matching conventions 
ttp = — 3^/^ai) while the second order term with its peculiar numerical coefficient is new. It is 
also straightforward to perform the strong coupling limit ai — )■ cxd for E{kf) in cq.(31). One 
finds E{kfY°°^ = (3A;j/10M)(— 1.948), corresponding to almost minus twice the free Fermi gas 
energy and therefore an instabihty of the system. The basic reason for this overly strong p-wave 
attraction is the enhancement factor 3 from the coherent sum of Hartree and Fock terms. In the 
case of a strong s-wavc interaction the number multiplying the free Fermi gas energy 3k^/10M 
is ^„ - 1 = -0.493 ~ -1.948 • 1/4. 



3.1 Particle-particle ladders only 

In this subsection wc study for comparison the resummation of only particle-particle ladders 
with a (spin-independent) p-wave contact-interaction. The factorization of multi-loop diagrams 
and the separation into transversal and longitudinal parts hold in the same way, while the 
resummation to all orders is accomplished now by a geometrical series. The resulting expression 
for the energy per particle reads: 



E{kf) 



72^ 
M 



djKK 



2/, 



+ 



n{aikf)-^ + F^{s, k) n{aikf)-^ + fI{s, k) 



(33) 



supplemented by a principal-value prescription for the double-integral. The transversal and 
longitudinal particle-particle bubble functions F^'^^{s,k,) occurring in the denominators have 
the form: 



8s2 



1 

3 + 

--(l-«;2) 
8s^ ' 



2\2 



o, 1 + S + K 1 — s"^ — 

In h 



X 



1 + S — K 



In 



16s3 

(1 + Sf - 
1 — — 



(34) 
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Figure 5: Energy per particle E{kf) from p-wave ladder diagrams divided by the free Fermi 
gas energy Sk^/IOM as a function of the dimensionless parameter aikf. 



FL{s,k) 



4 
3 



(1 - 



4s2 

1 + S + K 



(1 + sy - 



+ 



8s3 



1 



(35) 



These functions satisfy the relations -Pp^'"(s, + ReFp^''l(— s, k) = R±^\\{s, k) inside the region 
+ fi;^ < 1 of interest. As it must be the low-density expansion of E{kf) in eq.(33) agrees up 
to next-to- leading order with eq.(32). The additional hole- hole and mixed particle- hole ladders 
included in eq.(31) are of cubic and higher order in the scattering volume af. The strong cou- 
pling limit ai oo is also readily performed. One gets E{kf)^°°'^ = (3/cj/10M)(— 2.405), which 
amounts to about 23% more attraction than in the case of the complete ladder series. It is a 
general feature that hole-hole and mixed particle-hole ladders generate some moderate repul- 
sion. Again the enhancement factor 3 in comparison to the case of a strong s-wave interaction 
is recognizable here, C^p^ - 1 ~ -0.763 ~ -2.405 • 1/3. 

It is interesting to explore the dependences on density and coupling strength of the novel 
p-wave resummation results eqs. (31,33). For this purpose we plot in Fig. 5 the ratio of the 
interaction energy E{kf) to the free Fermi gas energy Skj/IOM as a function of the dimen- 
sionless parameter aikf. The full and dashed curve correspond to the complete ladder series 
and partial pp-ladder series, respectively. On the attractive side aikf > one observes that 
both curves stay closely together up to aikf ~ 1.5 and from there upwards the pp-ladder series 
develops somewhat more attraction. The behavior on the repulsive side aikf < is character- 
ized by a sharp peak of the full curve at aikf ~ —1.1 and an even higher peak of the dashed 
curve at aikf ~ —1.3. Beyond that a non-perturbative behavior sets in, such that both curves 
pass through zero at aikf ~ —1.35 and attraction gets generated from the repulsive scattering 
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volume al < 0. In order to join smoothly with the attractive side at infinity aikf = ±00, this 
attraction is more strongly pronounced for the pp-ladder series. 

As a side remark we note that s-wave and p-wave contact-interactions do not interfere in 
the ladder diagrams. The usual argument of rotational invariance (in momentum-space) does 
not apply in the medium. An sp-interference term introduces into the in-medium loop a factor 
that is odd under parity /—;■—/, but the integration regions defined by two overlapping or 
intersecting Fermi spheres in eqs.(5,6) are even under this transformation, and therefore it 
integrates always to zero. 

Let us add a final remark on the physical applicability of the resummation results derived 
in this work. In order to obtain the factorization of multi-loop diagrams (in a generalized 
sense) the input interactions have to be contact-interactions. Their iteration to all orders in 
the vacuum leads to a unitary S-matrix with phase-shifts given by eqs.(4,26). On the other 
hand the Wigner bound [21] due causafity states that a scattering phase-shift must not drop 
too fast with momentum, d5i{q)/dq > —R, where R measures the range of the interaction. 
In the case of (unresolved) contact-interactions the range ^ ~ vanishes practically. This 
consideration could provide an argument for the condition that the resummation results are 
physically applicable only for positive scattering length a > 0, positive scattering volume ai > 0, 
and negative effective range < 0. 

Appendix: Resummation with coupled channels 

In this appendix we supply the result for the complete resummation of in-medium ladder di- 
agrams due to an s-wave contact-interaction in a binary (i.e. two- component) many-fermion 
system. The prototype of such a system is nuclear matter made up from strongly interacting 
protons and neutrons. Low-energy nucleon-nucleon scattering is characterized by two different 
scattering lengths, the spin-singlet (/ = 1) scattering length Og ~ 19 fm and the spin-triplet 
(/ = 0) scattering length at ^ — 5.4fm, where / = 0, 1 denotes the total isospin. The many- 
body dynamics is richer and offers the possibility to introduce an isospin-asymmetry through 
different (variable) proton and neutron densities (or different Fermi momenta kp and A;„). Re- 
turning to the in-medium loop in Fig. 2 one sees that one has to deal now with four (partly 
coupled) states: {pp,nn,pn,np). The last two channels get decoupled through a 45° rota- 
tion to the isospin-eigenstates {pn ± np) / \f2. In this isospin-basis the momentum-independent 
contact-interaction (proportional to the scattering lengths a^^t) is diagonal and therefore the 
multi-loop diagrams factorize. As a new clement one encounters mixed proton-neutron loops 
with integration regions defined by eqs.(5,6), where one Fermi momentum kf is k^ and the 
other is The real part of this mixed in- medium loop is given simply by the sum of two 
independent proton and neutron terms, but a more complicated overlap of two Fermi-spheres 
|P ± / 1 < kpji and the surface \l \ = \q\ has to be evaluated for the imaginary part (see 
eq.(38) below). The resummation of the particle-hole ladder diagrams to all orders through 
an arctangent-function is also easily carried out, since for a 4 x 4 diagonal matrix one has the 
identity: arctan(diag(di, ^2, ^^3, ^^4)) = diag(arctan di, arctan ^2, arctan ^3, arctan ^4). In order 
to compute the energy density one rotates back to the particle basis and adds the Hartree and 
Fock contributions with their spin-factors 4 and —2. Putting all pieces together the final result 
for the energy per particle reads: 
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Figure 6: Support of the function $(P, g, kp, kn) in the (P, g)-plane pieced together by a triangle, 
a rectangle and a circular wedge. The rectangle has side-lengths \/2kp and — kp)/\/2. 



+ k„ arctan 



arctan 



lis, k) 



'flsK) ^ + TT ip(s, k) 
<^{P,q,kp, kn) 



{kp+k„)/2 



+ {2Tl)'\kpRp + knRn 



+ 3 arctan 



dPP'^ J dqq 

ain 



a^-i + (27r)-i(/i;pPp + A;„P„) 



(36) 



with the abbreviations Rp = R{P/kp, q/kp) and P„ = R{P/kn, q/kn), where R{s, k) and /(s, k) 
are given in eqs. (16,17). The integration boundaries are: 



gmin = max(0, P - kp) , g^ax = mm(^kp + P, y'(A;2 + k^)/2 - P^ 



(37) 



assuming for definiteness kn > /^p (i.e. an excess of neutrons over protons). The total nucleon 
density is p = Pp + pn = {kp + /c^)/37r^. The imaginary part function $(P, q, kp, kn) appearing 
in the numerator of the last two arctan-terms is defined piecewise by: 



$ = g , for P + g < fc, 



p ' 



*=4P 

4P 



kl-{P - qf] , hi kp < P + q < kn , \P - q\ < kj 



p ' 



kl + kl- 2(P2 + g2)l , for P + g > A;„ , P2 + g2 < (A;2 + A;2)/2 . (3^ 



As illustrated in Fig. 6 the support of this function consists of a triangle, a rectangle and a 
circular wedge pieced together. Note that P^g $(P, g, kp, kn) plays also the role of a weighting 
function (symmetric under P -H- g) in reducing the integral over two Fermi spheres |pi^2| < ^p,n- 
The same expression for the ^-function as in eq.(38) has been derived in appendix C of ref. [14]. 
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spin-asymmetry r| 

Figure 7: The (normal) Bertsch parameter in{v) of the unitary Fermi gas as a function of the 
spin- asymmetry parameter rj. 



The numerical implications of eq.(36) for the nuclear matter equation of state are discussed 
elsewhere [22]. The inclusion of corrections from the singlet and triplet effective ranges should 
also be obvious. The inverse scattering lengths a~l in the denominators of the arctan-tcrms 
in eq.(36) get then replaced by appropriate reciprocal O- functions. Note also the factor 3 in 
front of the last arctan-term, it ensures that the result E{kf,kf) for the isospin-symmetric 
system, kp — k^ — kf, becomes invariant under the interchange (a^ •<->■ at) of both scattering 
lengths. As a byproduct of eq.(36) wc consider the isospin- asymmetry energy A{kf). It is 
obtained by setting the Fermi momenta to kn^p = A:/(l ± 5)^^^ and expanding to quadratic 
order, E{kp, kn) = E{kf) + S'^A{kf) + 0{6^), in the isospin-asymmetry parameter 5. One finds 
up to second order in the scattering lengths Gs/. 

k"^ ( k S/c^ 1 

^^^^^ = 6M + i^^""' ~ + 5^["'^^ a?(2 + ln2)] + . . . | , (39) 

where the first term k'j/GM arising from the kinetic energy has been included for orientation. 
This result is implicitly contained in ref.[23] via the parameterizations Gg^t = g\Mli,o/{^T^fl) 
of the scattering lengths. As a special case of an asymmetric system let us finally consider 
the one-component (unitary) Fermi gas in the spin-unbalanced configuration. With one singlet 
scattering length a there is interaction taking place only in the mixed spin-state (ti — it) / 
and the energy per particle reads: 

ymax T / I-, 1 1 \ 

^min 

with ^max, Q'min, -Rf' ^"^^ Q') ^t? ^4-) defined in complete analogy to the proton-neutron 
system before. Note that the density of the spin-unbalanced Fermi gas is p = (A;| -|- /c|)/67r^. 
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We can use eq.(40) to determine, in the unitary limit a — )■ oo, the dependence of the (normal) 
Bertsch parameter ^niv) the spin- asymmetry r] of the spin-unbalanced system. Setting the 
Fermi momenta of the spin- up and spin-down fermions to A;^,^, = ± r^)^/^, the quantity 
iniv) — 1 is given by the ratio of the interaction energy E{k^, k\) to the kinetic energy 3A;j[(l -|- 
jy)5/3 _|_ — ri)^l^]/2QM . The full line in Fig. 7 shows the (normal) Bertsch parameter inij]) 
as a function of the spin-asymmetry parameter rj. The result is quite close to the parabolic 
approximation ^nij]) = (1 + '7^)/2 connecting the two boundary points {r],^n) = (0, 1/2) and 
(r?,a) = (l,l)- 
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